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Abstract The model of current-carrying string loop oscillations is tested to 
explain the special set of frequencies related to the high-frequency quasiperi¬ 
odic oscillations (HF QPOs) observed recently in the low-mass X-ray binary 
XTE J1701-407 containing a neutron star. The external geometry of the neu¬ 
tron star is approximated by the Kerr geometry, introducing errors not ex¬ 
ceeding 10 % for slowly rotating massive neutron stars. The frequencies of the 
radial and vertical string loop oscillations are then governed by the mass M 
and dimensionless spin a of the neutron star, and by the dimensionless param¬ 
eter lo describing combined effects of the string loop tension and its angular 
momentum. It is explicitly demonstrated that the string-loop oscillation model 
can explain the observed kHz frequencies for the neutron star parameters re¬ 
stricted to the intervals 0.2 < a < 0.4 and 2.1 < M/ M 0 < 2.5. However, the 
stringy parameter uj cannot be the same for all the three HF QPO observa¬ 
tions in the XTE J1701-407 source; the limits on the acceptable values of oj 
are given in dependence on the spacetime parameters M and a. 

Keywords XTE J1701-407 • QPOs • string-loop 
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1 Introduction 

The axisymmetric current-carrying string loops are governed by their tension 
and angular momentum. Tension of the string loops prevents their expansion 
beyond some radius, while their worldslieet current introduces an angular mo¬ 
mentum preventing them from collapse. Originally, the cosmic strings were 
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introduced as remnants of the phase transitions in the very early universe 
[84], or strings represented as superconducting vortices were considered in 
[88]. However, the current-carrying string loops could represent also plasma 
exhibiting a string-like behaviour due to dynamics of the magnetic field lines 
[53,20], or due to the thin flux tubes of magnetized plasma simply described 
as ID strings [54,21,22,40]. 

In recent studies it has been demonstrated that the current-carrying string 
loops moving axisymmetrically along the symmetry axis of the Kerr or Schwarz- 
schild-de Sitter black holes can have interesting direct astrophysical appli¬ 
cations [33,36,62,37]. For example, the transmutation effect related to the 
current-carrying string loops, i.e., transmission of their oscillatory internal en¬ 
ergy into energy of the translational motion, causes an outward-directed accel¬ 
eration of the string loops in the gravitational field of stars or compact objects, 
as neutron stars, black holes, or naked singularities [33,62,63,37]. Such an ef¬ 
fect can be important also for the electrically charged string loops moving in 
combined external gravitational and electromagnetic fields [80,81]. 

Since acceleration of the string loops in the deep gravitational field of 
black holes or naked singularities can be extremely efficient, leading to ultra- 
relativistic escaping velocities of the string loops [63,37], the string loop trans¬ 
mutation can serve as a model of the formation of the ultra-relativistic jets 
observed in active galactic nuclei or Galactic microquasars. The string loop 
transmutation model of collimated relativistic jets can thus be considered as 
an alternative to the standard model based on the Blandford-Znajek process 
[10], or recently proposed model of geodesic collimation of ultra-relativistic 
particles [24,25]. For acceleration of jets escaping from the active galactic nu¬ 
clei the cosmic repulsion can be also relevant behind the so called static radius 
[62,60,61]. 

Quite recently, it has been demonstrated that the string loops can be as- 
trophysically relevant in situation inverse to the transmutation effect. Small 
oscillations of the string loops around stable equilibrium positions in the equa¬ 
torial plane of the Kerr geometry can be considered in the lowest approxima¬ 
tion as two uncoupled linear harmonic oscillators governing the radial and 
vertical oscillations of the string loop - the higher-order terms determine non¬ 
linear phenomena and subsequent transition to the quasi-periodic and chaotic 
oscillatory motion [37]. The frequencies of the radial and vertical harmonic 
oscillations of the string loops are relevant also in the quasi-periodic stages of 
the motion and their radial profiles were given and discussed in [64], It has 
been shown that frequencies of the string loop harmonic or quasi-harmonic 
oscillations can fit twin HF QPOs observed with frequency ratio 3 : 2 in three 
Galactic microquasars GRS 1915+105, XTE 1550-564, GRO 1655-40, i.e, low- 
mass X-ray binary (LMXB) systems containing a black hole [64], 

The radial profiles of the string loop oscillations qualitatively differ from 
those related to the radial and vertical oscillations of the geodesic, test particle 
motion in the Kerr geometry, especially there is a crossing point of the radial 
and vertical frequencies in the Kerr black hole spacetimes for the string loop 
oscillations, while for the test particle oscillations such a crossing is possible 
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only in the Kerr naked singularity spacetimes [78,70]. This property of the 
radial profiles of the radial and vertical frequencies of the string loop oscilla¬ 
tions is applied here to test the string loop oscillation model [64] for fitting the 
special frequency set of HF QPOs observed recently in the source XTE J1701- 
407 that is a LMXB system containing a neutron star [50] . For simplicity, we 
assume that the external geometry of the neutron star can be appropriately 
approximated by the Kerr geometry. 


2 HF QPOs observed in the XTE J1701-407 source 

The LMXBs can contain a black hole or a neutron star. The neutron star 
LMXBs are classified as Z sources if they demonstrate high luminosities (0.5 — 
1.0 of Eddington luminosity TEdd), and atoll sources if they demonstrate low 
luminosities (0.01 — 0.5 Lecm) [29,83]. 

The source XTE J1701-407 has been discovered by the Rossi X-ray Timing 
Explorer (RXTE) at June 2008 [45]. Since thermonuclear X-ray bursts were 
observed at the XTE J1701-407 source [19], we can conclude that it contains 
an accreting neutron star. The discovery of HF QPOs at the XTE J1701-407 
source was first reported in [59] where it was mentioned that it is one of the 
least luminous atoll sources, having the X-ray luminosity Lx ~ O.OILecm- 
Later, it has been discovered that the source is a peculiar neutron star LMXB 
displaying all various types of the neutron star LMXB subclasses, i.e., Z-class, 
bright atoll and atoll [30] 

Recently, a detailed analysis of the HF QPOs in the source XTE J1701- 
407 has been reported in [50] , being based on the whole ansamble of 58 RXTE 
observations. The results are rather unexpected and very interesting, since 
a very special set of frequencies has been discovered at three of the whole 
ansamble of the RXTE observations. In one of the three observational events 
a single HF QPO has been detected at a characteristic frequency [50] 

/(A)l = /(A)u = H53±5Hz. (1) 

In the other two observations twin HF QPOs has been detected at character¬ 
istic frequencies [50] 

/( B)L = 740 ± 5 Hz, /( B)U = 1112 ± 17 Hz, (2) 

/( C )L = 740 ± 11 Hz, / (C)U = 1098 ± 5 Hz, (3) 

where we use the index U for the upper and the index L for the lower of the 
twin frequencies observed simultaneously. The other details of the detected HF 
QPOs, namely the fractional rms amplitudes and the quality factors, are pre¬ 
sented in [50] . We do not discuss them here, since we concentrate our attention 
on the observed frequencies only. 

During the observations when the one-peak HF QPO has been detected, a 
low-frequency QPO has been also detected at fi ow = 30.4 ±0.3 Hz [50] that is 
not considered in the following as we restrict our attention to the HF QPOs. 
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However, we plan to study the possibility that the low-frequency QPO could 
be related to ’’rocking” of a string loop around the symmetry axis. 

Note that the HF (kHz) QPOs do not always appear as twin peaks, but 
only one peak is sometimes observed in the neutron star systems. If only one 
kHz peak is detected, it could mean that the other peak is not detected for 
physical reasons, or because of limitations of detectors, such as low signal-to- 
noise. The physical reasons could block production of the second peak, or, as 
we assume here, frequencies of the two oscillatory modes coincide. 

The frequency set corresponding to the HF QPOs detected in the XTE J1701- 
407 source is really of a special character among the HF QPOs observed in the 
atoll sources, since a single HF QPO is combined with two twin HF QPOs hav¬ 
ing the frequency ratio 3 : 2 typical for microquasars. We shall use the special 
character of the radial profiles of the frequencies of the harmonic radial and 
vertical oscillations of axially symmetric string loops, namely the coincidence 
of these frequencies at a special radius [64], in order to fit quite naturally the 
single HF QPO observed in the source XTE J1701-407. 


3 Neutron stars and Kerr geometry 

In our test of the string loop model for the case of the neutron star system 
XTE 1701-407 we apply the Kerr geometry to describe the external geometry 
of the neutron star, as the radial profiles of the radial and vertical string loop 
oscillation frequencies are given by analytic formulae in this geometry and 
enable relatively simple and elegant fitting procedure. However, we have to 
clear up the conditions enabling using of such a strong simplification for the 
external geometry of the neutron star. 


3.1 Hartle-Thorne theory of rotating neutron stars 

The rotating neutron stars can be conveniently described by the Hartle-Thorne 
geometry that is governed in its exterior part by three parameters: mass M, 
dimensionless spin a, and dimensionless quadrupole moment q. The Hartle- 
Thorne theory can be applied for rotating neutron stars with rotation fre¬ 
quency smaller than the mass-shedding frequency (/mass shedding ~ 1100 Hz) 
and dimensionless spin a < 0.4 [27,28,18,47,82]. For a < 0.4 and q/a 2 = 1, 
the external Hartle-Thorne geometry [28] coincides with the Kerr geometry 
with precision on the level of ~ 1 %. For 1 < q/a 2 < 3 these two geome¬ 
tries have very similar radial profiles of the metric coefficients, giving thus for 
astrophysical phenomena occurring in vicinity of such rotating neutron stars 
predictions with errors on the level lower than ~ 5 % [74,75,76,7,9,79]. 

Recently a very important result has been obtained in [82] where it has 
been demonstrated that for the whole variety of considered realistic equations 
of state enabling existence of massive neutron stars with M > 2M 0 , the Hartle- 
Thorne theory predicts q/a 2 < 2, if the mass of the neutron star with a given 
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rotation frequency is close to the maximum allowed for a given equation of 
state. We can see in Figure 3 of [82] that near the mass-maximum given by 
equations of state allowing for M max > 2M 0 , the value of the characteristic 
parameter can be as low as q/a 2 ~ 1.5. Therefore, such near-maximum-mass 
neutron stars can be properly described by the Kerr black hole geometry. The 
precision of the agreement of the external Hartle-Thorne and Kerr geometries 
is sufficiently high (with relative difference lower than 3 %), if the conditions 
a < 0.4 and q/a 2 < 2 are satisfied [75,82]. 

For masses substantially lower than the maximum allowed by the realis¬ 
tic equations of state, the parameter q/a 2 can be relatively high - for M ~ 
0.7 M max there is q/a 2 —>• (5 — 10) in dependence on a concrete equation of 
state [82]. We can see in Figure 3 of [82] that for the so called ’’canonical” 
neutron star mass M = 1.4Mq, there is q/a 2 —> (3 — 7) for all of the con¬ 
sidered realistic equations of state. In such cases, dependence of astrophysical 
phenomena on the neutron star parameters M, a, q can be very complex as 
demonstrated in [79] for the case of the position of the innermost stable cir¬ 
cular geodesic (risco) in the Hartle-Thorne spacetime. The difference of the 
nsco in the Hartle-Thorne and Kerr spacetimes having identical parameters 
M and a reaches tens of percent and the dependence of risco on the spacetime 
parameters differs even qualitatively for any value of spin a < 0.4 and for large 
values of the parameter q giving q/a 2 > 3. On the other hand, the relative 
differences are less than 10% for spin a < 0.4 and q/a 2 < 3 - see Figure 1 in 
[79]. 

We can summarize that for the near-maximum mass neutron stars the sit¬ 
uation is common and relatively simple for all equations of state if the spin 
parameter a < 0.4, since the quadrupole moment parameter satisfies the con¬ 
dition q/a 2 < 2 implying similarity of the Hartle-Thorne and Kerr geometries 
and similar character of astrophysical phenomena in these spacetimes, with rel¬ 
ative differences smaller than 5 %. The Hartle-Thorne metric is quadratic in 
the dimensionless spin a, and for a ~ 0.4 it introduces truncation errors of the 
order 0(a 3 ) ~ 6 % implying total errors introduced by assumption of the Kerr 
background geometry about 10 %. Therefore, more detailed and accurate anal¬ 
ysis of the oscillatory motion, based on the Hartle-Thorne model or numerical 
codes of fast-rotating neutron stars (RNS [31,38] or LORENE [32,12] codes) 
will be necessary in future studies. For masses significantly smaller than the 
maximal mass related to a given equation of state, the Hartle-Thorne geome¬ 
try depends strongly on the equation of state having the parameter q/a 2 > 3; 
in such a case it differs significantly (qualitatively) from the Kerr geometry. 
Approximation of a neutron star exterior by the Kerr geometry is irrelevant 
in such situations. 


3.2 Theoretical restrictions on the mass and spin of neutron stars 

Predictions of the string-loop oscillation model for the XTE J1701-407 source 
containing a neutron star have to be confronted with the theoretical con- 
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straints on mass and spin of neutron stars. The upper theoretical constraint 
on the neutron star mass based on realistic equations of state of the subnuclear 
matter reads 

M < M max NS = 2.5 M 0 (4) 

[1,2,17]. Note that the limit governed by the causality condition gives M maxNS ~ 
3Mq [17], and the equation of state based on the relativistic mean-field theory 
[11] allows for M max NS ~ 2.8M 0 . 

The upper theoretical limit on the neutron star spin, based on the fully 
general relativistic set of numerical codes for numerical relativity, LORENE 
[32], reads a max NS ~ 0.7, as demonstrated in [44], and corresponds to the 
maximal rotational frequency of the neutron stars, / rot = /mass shedding- For 
/rot ~ 0.6/ masa shedding corresponding to the observed neutron stars with 
fastest rotation, the LORENE codes predict the neutron star spin a < 0.4 
- see Figure 3 in [44]. Surprisingly, the numerically obtained fully general rel¬ 
ativistic, fast rotating, ’’LORENE” external geometry, and the Hartle-Thorne 
external geometry obtained in the ’’slow rotation” approximation of the neu¬ 
tron star models, differ for the equations of state considered in [82] maximally 
on the level of 5%, if the near-maximal mass neutron stars with the dimen¬ 
sionless spin a < 0.4 are considered. Their difference is substantially higher, 
> 10%, for the spin a > 0.5. Therefore, we put the theoretical limit on appli¬ 
cability of the Hartle-Thorne model of neutron stars to be 

a <0.4. (5) 

For quark (strange) stars, the maximal mass constraint is smaller because 
of softer equations of state assumed in modelling the strange stars, but masses 
around M max q ~ 2M 0 are allowed [26,44]; in some special conditions the 
maximum mass can be substantially higher, namely M max q ~ 2.7M 0 [17]. 
On the other hand, slightly superspinning states of strange quark stars with 
OmaxQ > 1, exceeding the black hole limit a = 1, are allowed as reported in 
[44] . Such a difference between the maximal spin of neutron stars and strange 
stars can be explained by the strong nuclear force acting in binding the strange 
stars. Here we consider only the neutron star models. 

We use the theoretical restrictions (4) and (5) as limits on the mass and 
spin of the XTE J1701-407 neutron star, given by the string loop oscillation 
model applied to fit the frequencies of the HF QPOs in the three observational 
events reported in [50]. 


4 Oscillations of current-carrying string loops 

An axisymmetric current-carrying string loop in a given axially symmetric and 
stationary spacetime with metric g a p has been discussed in detail in [33,37]. 
Oscillations of such string loops have been studied in [64] . Here we give a short 


overview. 



String loop oscillation model in a neutron star source 


7 


4.1 Dynamics of axisymmetric string loops 


The worldsheet of a string loop is given by the spacetime coordinates X a (a a ) 
with a = 0,1,2, 3 that are functions of two worldsheet coordinates <r a with 
a = 0,1. The induced metric on the worldsheet reads h a b = g a pX ,“Xj^ where 
□ u = dU/da. The string loop current is governed by a scalar held <p(cr a ). 
Dynamics of the string loop is determined by an action with Lagrangian £ = 
— (g + h ab ip\ a ip\b)\/—h where <p }a = j a determines the current of the string 
and g > 0 reflects the string tension [84,88,33,36]. Axisymmetry of the string 
loop in accord with axisymmetry of the spacetime means that the scalar held 
tp = j a <7 + j T r, where j a and j T are constant components of the current and 
represent the angular momentum of the string loop. We dehne the angular 
momentum parameters by the relations J 2 = j 2 + j 2 and uj = — j a /j T [37]. 

Dynamics of the axisymmetric string loops is then given by the worldsheet 
stress-energy tensor expressed through the tension g > 0 and the parameters 
J > 0 and ui £ (—1,1) that are related to the current on the string loop 
and determine its angular momentum. For a given string loop, the tension 
and both the angular momentum parameters are constants of the motion. 1 
The worldsheet stress-energy tensor of the string loop X ab reads (see [64] for 
details) 


S TT 


J2 

-b g, X' 

94 ><!> 



2 ujJ 2 

g^(l + w 2 )' 


( 6 ) 


The string loop momentum components are given in the standard way as 
P'n = where dot means derivative with respect to the proper time r. 
Dynamics of the string loops is then governed by the Hamilton equations [42, 
43,23] with the Hamiltonian given in the form [42,37] 

H = \g at) P a P? + \g H [(X TT ) 2 - (P Ta ) 2 ] ■ (7) 


The equations of evolution of axisymmetric string loops are similar to the 
equations governing evolution of toroidal magnetic fields and plasmas in a 
given background [53,20,33,36]. Such composites of toroidal magnetic fields 
nad plasmas can be generated by the kinetic dynamo effect [21], The oscillating 
string loops can represent a simplification of oscillating toroidal structures of 
magnetized plasmas, reflecting by their tension the effect of trapped magnetic 
toroidal field, while their angular momentum reflects the role of the matter 
content of such toroidal configurations. The oscillating string loop model can 
be thus considered as complementary to models of oscillating perfect-fluid tori 
reflecting the role of pressure gradients [51,58] where the effect of tension is 
not taken into account. We expect that physical parameters of a thin flux tube 
of magnetized plasma can be approximated by the string loop parameters J, ui. 


1 For an alternate approach to the modelling of strings see the vorton model of [16,15], 
or strings with time-varying tension [87]. A detailed information and review can be found 
in [84], 
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4.2 Energy boundary function 


The stationarity and the axial symmetry assumed for the spacetime imply 
existence of two conserved quantities - the string-loop energy E and axial 
angular momentum L that are given by 

— E = Pt, L = = g < p ( pE aT = —2 j T j a . (8) 


The axisymmetric current-carrying string loops are thus characterized by the 
total angular momentum J and the axial angular momentum L. The case 
u> < 0 (u) > 0) represents string loops with negative (positive) axial angular 
momentum L < 0 (L > 0); u = 0 represents the string loops with L = 0 [33, 
37,64], 

The Hamiltonian can be written as a sum of the dynamic and potential 
parts 

H = H D + H P = ^g rr P 2 + lg ee P-f + H P (r, 9). (9) 

Using the energy E and the angular momentum parameters J and u>, the 
potential part of the Hamiltonian reads [64] 


H P 




1 9U ( E 9t<)> 2 J 2 u 

^(gttg^-gU) ^ 9 ** u2 +1 


2 


( 10 ) 


The boundary of the string loop motion is given by the condition Hp = 0 that 
implies the relation [37,64] 


E = E h {r, 9) = y/ g 2 t(j> - g tt g^ E TT - g t4 E aT . (11) 


The energy boundary function Eb(r,9) governs the dynamics of the string 
loops, serving as an effective potential of their motion. 

We express the energy boundary function in the Kerr metric using the 
motion constants J and w, making the rescaling E^/fi —> E b and J/y/JI —> J 
- such a choice of “units” is equivalent to setting the string tension g = 
1. The energy boundary function then, in the standard Boyer-Lindquist r, 9 
coordinates [14], takes the form [64] 

_ , „ T . 4awJ 2 r r— ( J 2 R 2 . \ 

E b (r, 8;a, J,u) = + ^ + 8 ‘” ( ’) ’ ( ’ 

with 

G(r,9;a)= (a 2 + r 2 ) R 2 + 2a 2 r sin 2 (0). (13) 

The Kerr metric is governed by the characteristic functions 

R 2 = r 2 + a 2 cos 2 0, A = r 2 -2Mr + a 2 , (14) 


where a denotes spin and M mass parameters of the Kerr spacetimes. In 
the following, we shall use for simplicity the dimensionless radial coordinate 
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r/M —» r, dimensionless time coordinate t/M —► t and dimensionless spin 
a/M —> a; this is equivalent to using of M = 1 in the metric tensor. We will 
return to the dimensional quantities in the Section 5. 

Detailed discussion of the properties of the energy boundary function Eb(r, 6) 
is presented in [37] for both the Kerr black hole and naked singularity space- 
times. Here we focus on its properties in the black hole spacetimes that can 
be relevant for rotating neutron stars as demonstrated in [82,74] - the local 
extrema of the energy boundary function can be then located in the equatorial 
plane only. 

The stationary points, i.e., the local extrema of the energy boundary func¬ 
tion Ei,(r\ a, J, w), governing the equilibrium positions of the string loops in 
the equatorial plane (6 = 7r/2), are determined by the function Jg(r;a, u>) 
defined by the relation [37,64] 


(r — 1) (w 2 + l) H 2 
4 ajj\/~A ( a 2 + 3r 2 ) + (w 2 + 1) F’ 


(15) 


where 


H(r; a) = a 2 (r + 2) + r 3 , F(r;a) = (r — 3)r 4 — 2a 4 + a 2 r (r 2 — 3r + 6) . (16) 

A detailed discussion of the properties of the energy boundary function 
Eb(r; a, J,lj) and the string loop motion can be found in [37,64]. We have to 
concentrate on the situation when for a string loop with fixed values of the 
angular momentum parameters J and ui closed E = const sections of the effec¬ 
tive potential (energy boundary function) occur around a stable equilibrium 
position of the string loop given by the equation 

J 2 = (17) 

The stable equilibrium position corresponds to the minimal energy related to 
the string loop with the angular momentum parameters J, u>. Around such sta¬ 
ble equilibrium positions, small oscillations of string loops occur, if their energy 
slightly exceeds the minimal value. The analysis of the oscillatory motion of 
string loops around their stable equilibrium positions, using the perturbative 
treatment of the Hamiltonian, can be found in [63,37,64]. 


4.3 Frequency of the string-loop radial and vertical oscillatory modes 


We consider string loop harmonic oscillations around a stable equilibrium po¬ 
sition with fixed coordinates ro and 9q = 7r/2 corresponding to the local min¬ 
ima of the energy boundary function. For the radial and latitudinal oscillatory 
string loop motion in the Kerr spacetimes, the angular frequencies related to 
distant observers are given by [64] 


f? 2 (r) = 




E(ex) 


2auV A (a 2 + 3r 2 ) + u> ( a 2 r 3 — a 2 A + r 5 — 2?- 4 ) 


2rH 2 


2 auj (a 2 + 3r 2 ) + wa/A (r 3 — a 2 ) 


, ( 18 ) 









10 


Z. Stuchlfk, M. Kolos 


a (2 lo\]A — ati'j [a 2 (2 — 3r) — 3r 3 ] + r 2 a> [a 2 (r — 6) + r 3 ] 


ttg(r) = 


r 2 HA -V2 


2 aco (a 2 + 3r 2 ) + w-\/z\ (r 3 — a 2 ) 


(19) 


where ui = (w 2 + l) and 

J E (ex)( r ) = i?w(r — 1) (6a 2 r — 3a 2 r 2 — 6a 2 — 5 r 4 + 12r 3 ) 

+w [2F(a 2 + 3r 2 )(l - r) - FH] + 8auy/A(r - l)(o 2 + 3r 2 ) 2 
+4clujA~ 1 / 2 H [(a 2 + 3r 2 ) (A — (r — l) 2 ) — 6 rA(r — 1)] . (20) 

The function JE(ex)( r > a > w ) governs the local extrema of the function Je{t", a, u>). 
Its zero points determine the marginally stable equilibrium positions of the 
string loops - at the zero points the frequency of the radial oscillatory modes 
of the string loops vanishes. 2 The string-loop oscillations are possible only if 
the stable equilibrium positions of the string loops are allowed. The conditions 

■^E(ex) = 0 and J E > 0, (21) 

satisfied simultaneously, put the limit on validity of the formulae giving the 
angular frequencies of the radial and vertical oscillations - for details see [64], 
The radial profiles of the frequency of the radial and vertical modes of the 
string loop harmonic (or quasi-harmonic) oscillations has been discussed in 
[64]. A very complex situation occurs for the case of near-extreme Kerr black 
hole spacetimes with a > 0.99, similarly with the situation of the test particle 
and fluid motion related to the so called Aschenbacli effect [3,72,71]. In the 
Kerr spacetimes with dimensionless spin a < 0.7 which has been shown to be 
the upper limit on the spin of neutron stars, independently of the equation of 
state [44] , the properties of the radial profiles of the string loop oscillations are 
much simpler, as demonstrated in Figure 1 for four characteristic values of the 
Kerr spin parameter a = 0,0.2, 0.4,0.7; in the Schwarzschild spacetime (a = 0), 
the situation is degenerated, since both the frequencies are independent of the 
parameter u>. The range of the radii where the oscillating string loops can occur 
due to possible variance of the stringy parameter u> increases with increasing 
spacetime dimensionless spin parameter a, as illustrated in Figure 2. 


5 String-loop oscillations as a model of HF QPOs in the XTE 
J1701-407 source 

We realize fitting of the frequencies of the kHz QPOs observed in the neutron 
star source XTE 1701-407 using the analytical radial and vertical frequency 
relations of the string loop harmonic oscillations in the Kerr geometry. Our 
preliminary results [65] demonstrate that the relative difference of the fre¬ 
quency relations in the Kerr geometry with a < 0.4 and the Hartle-Thorne 

2 Exactly the same situation occurs in the case of the radial epicyclic motion of test 

particles [78]. 
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Fig. 1 String-loop oscillatory frequencies u r (thin curves) and uq (thick curves), calculated 
for the Kerr metrics with M = 2.5Mq. Their radial profiles are illustrated for values of 
dimensionless spin a = 0, 0.2, 0.4, 0.7 that are characteristic of our study of neutron star 
system. We demonstrate extension of the frequency radial profiles for the complete range of 
the string loop parameter lu E (—1,1). The vertical frequency curves are restricted to the 
region of existence (zero point) of the corresponding radial frequency curves - the relevant 
region is greyed. 




/(A) 

/(B) 

/(C) 

case 1 

v Y : u e 

1:1 

3:2 

2:3 

case 2 

v x : u e 

1:1 

2:3 

3:2 

case 3 

u r : u e 

1:1 

3:2 

3:2 

case 4 

u r : u e 

1:1 

2:3 

2:3 


Table 1 Four possible combinations of HF QPOs observed in the XTE J1701-407 source. 


geometry with a < 0.4 and q/a 2 < 2 do not exceed 5 %, giving thus sufficient 
accuracy to approve the application of the analytic Kerr formulae instead of 
complex, numerically treated Hartle-Thorne formulae. We also assume rele¬ 
vance of the resonant phenomena of the string loop oscillations at the radii 
where 3 : 2 and 1 : 1 frequency ratios occur. 
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Fig. 2 The resonant radii where the radial and vertical frequencies of the string-loop 
oscillations have rational ratio 3 : 2, 1 : 1 and 2 : 3, in dependence on the Kerr spin 
parameter a and the complete range of the string-loop parameter lj\ the case of uj = 0 is 
given by thin lines, the upper thick lines are given for uj = — 1 , the lower thick lines for 
uj = +1- The loci of the outer horizon (thick) and radii r = 3,4,5 (dashed) that represent 
possible radii of the neutron stars are also presented. Clearly, the string loop oscillations are 
allowed only above the surface of the neutron stars. 


5.1 Fitting observational frequencies by frequencies of the oscillatory modes 
of string loops 


The frequencies of the string loop oscillations measured by the distant ob¬ 
servers are in the dimensional form given by 


v {t,8) 


1 c 3 
2tt GM 




( 22 ) 


The frequencies of the harmonic string loop oscillations in the Kerr geometry 
have been compared to preliminary results of the analysis of the string loop 
oscillations in the Hartle-Thorne spacetimes [65] in order to justify the applica¬ 
tion of the Kerr geometry. The comparison has been done for the cases where 
we expect small differences in the frequency radial profiles, namely for spin 
a = 0.4 and the parameter q/a 2 = 1 and q/a 2 = 2. For comparison the case 
of q/a 2 = 10 is included. We can see in Figure 3 that the relative difference is 
really low for the cases q/a 2 = 1 and q/a 2 = 2, and high for q/a 2 = 10. For 
the neutron stars with near-maximal mass, dimensionless spin a < 0.4, and 
low parameter q/a 2 < 2, the Kerr formulae for the harmonic frequencies vg 
and t'r imply relative difference to the corresponding Hartle-Thorne formulae, 
I^Kerr — ^ht|/^ht, lower than 5 %. 

In the LMXB systems containing a black hole or a neutron star, some 
of the HF QPOs come in pairs of the upper and lower frequencies (iqj, vq/) 
of twin peaks in the Fourier power spectra that are sometimes combined with 
individually detected HF QPOs. In the LMXB systems containing a black hole, 
so called microquasars, the twin HF QPOs occur at fixed frequencies that are 
usually at the exact 3 : 2 ratio [46] . On the other hand, in the LMXB systems 
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Fig. 3 Radial profiles of the string loop harmonic oscillation frequencies vg and u r , related 
to the distant observers, are given for the Hartle-Thorne and Kerr spacetimes, and for the 
string loop parameter u £ {—1,0,1}. Both the Hartle-Thorne and Kerr spacetimes are 
rotating with the dimensionless spin a = 0.4, while the quadrupole moment q of the Hartle- 
Thorne spacetime takes the values q £ {0.16,0.32,1.6}. For Hartle-Thorne spacetimes, the 
string loop oscillation frequencies are given only numerically. In the case of q = 0.32 ( q/a 2 = 
2), the relative difference I^Kerr — ^htI/^ht of the Hartle-Thorne and the Kerr frequency 
radial profiles taken at the radius where the ratio vg : v T = 1 : 1, is not exceeding 5 %. 


containing a neutron star, large scatter of the twin HF QPOs is observed, and 
their frequency ratio varies across the frequency ratio 3 : 2 down to 5 : 4 ratio 
in the atoll sources [6,8,74,85,86], while it varies in the range 4 : 1 to 2 : 1 in 
the Z-source Circinus XI [13,75]. 

Usually, the Keplerian orbital and epicyclic (radial and latitudinal) fre¬ 
quencies of geodetical circular motion [56,78,39,66,52,70,49,5] are assumed 
in models explaining the HF QPOs in both black hole and neutron star sys¬ 
tems. Alternatively, the oscillations of tori [51] or tilted oscillating discs [34] 
are considered as explanation of the HF QPOs. The frequencies of the oscil- 
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lating tori and tilted discs are related to the orbital and epicyclic geodesic 
frequencies, if the oscillations are governed mainly by the gravity of the black 
hole or the neutron star [68]. 

The assumption that resonances of oscillatory modes of Keplerian accretion 
discs have to be relevant seems to be confirmed due to the properties of the HF 
QPOs observed in the three Galactic microquasars GRS 1915+105, XTE 1550- 
564, GRO 1655-40, where the 3 : 2 frequency ratio has been observed [35,73]. 
However, the observed frequencies of the HF QPO in the three microquasars 
cannot be explained by a unique model based on the frequencies of the geodesic 
epicyclic motion [77,4]. In the LMXBs containing neutron stars the fitting of 
the twin HF QPOs is realized by the relativistic precession model [55,57] or 
its variations [75]. However, the quality of the fits is not satisfactory for the 
models based on the geodesic orbital and epicyclic frequencies [76], although 
in some sources (4U 1636-53) the Resonant Switch model [67,68] can improve 
the fits significantly without introducing a new free parameter in addition to 
the spacetime parameters [69]. 


The quasi-periodic character of the oscillatory motion of current-carrying 
string loops trapped in a toroidal space around the equatorial plane of the 
central strong gravity object suggests application to the HF QPOs observed 
in the LMXB systems because of the same mass scaling of the frequencies 
of the radial and vertical oscillatory modes of string loops as those of the 
epicyclic geodesic frequencies, and comparable magnitudes of the frequencies 
of the string loop and epicyclic oscillations [64]. The string loop oscillation 
models has been successfully applied for the black hole systems in [64], here 
we test the model for the special kHz frequency set observed in the neutron 
star system XTE 1701-407. 

For a given twin HF QPOs observed in a given source, we have to consider 
fixed values of the string parameter to and the spacetime parameters M and 
a. If several twin HF QPOs are observed in the source, the spin and mass 
parameters have to be fixed, but the string loop parameter w can be varied, 
as different twin frequency observations could be generated by different string 
loops that could be created and decayed successively with different parameter 
u) reflecting locally different conditions in the accretion discs of the source. 
Therefore, the string-loop oscillation model naturally introduces a possibility 
of significant scatter in distribution of frequencies of the twin HF QPOs. Such 
a scatter can occur because of shifting of the created string loops along the 
accretion disc, or due to varied conditions at a fixed radius of the disc. On 
the other hand, the resonance phenomena can be relevant for the behaviour 
of the oscillating string loops as indicated by the Kolmogorov-Arnold-Moser 
(KAM) theory [48,64], or some resonance phenomena could be relevant even 
for creation of the string loops, selecting thus some special radii related to the 
resonant phenomena. 
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case 3 



Fig. 4 Restrictions on the mass M and spin a parameters of the neutron star in the 
XTE J1701-407 source implied by the string loop oscillation model applied to the three 
observational events of HF QPOs at the source. We assume that the three observational 
events occur at the resonant points of the radial and vertical string loop oscillations. Four 
different cases of the combinations of the resonant points related to the three observational 
events are possible - see Table 1. In the fitting procedure related to the upper observed 
frequencies, we have included the mean values of the observed frequencies (thick) and both 
the lower and upper values given by measurement errors (thin), see (1-3). The upper branches 
are for parameter co = +1 and the lower branches for parameter to = — 1, allowed regions 
of the spacetime parameters M, a are hatched. The most promising is the case 4, where 
we consider the 1 : 1 resonance with frequency = 1153 Hz combined with two 2 : 3 
resonances with frequencies V 2-.3 = 1088,1112 Hz. 


5.2 Fitting the kHz frequencies observed in XTE J1701-407 


The special set of the HF QPOs observed in the XTE J1701-407 source in¬ 
dicates application of the string loop oscillation model assuming a resonance 
entering the play as in the microquasars. The presence of the single HF QPO 
observed in the XTE J1701-407 neutron star system can be well attributed to 
the fundamental property of the radial profiles of the string loop oscillatory 
radial and vertical modes, giving for each value of the stringy parameter to a 
resonant position with frequency ratio u T : = 1:1. 
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case 1: 



a = 0.350 

a = 0.5 

a = 0.7 

M/M 0 

2.47 

2.38 — 2.61 

2.28 — 2.85 

^3:2 

-1 

-1 - -0.27 

-1 — 0.03 

^2:3 

-0.02 

-0.12 — 0.19 

-0.21 — 0.37 

Wl;l 

1 

0.30 — 1 

0.11 — 1 

case 2: 


a = 0.389 

a = 0.5 

a = 0.7 

M/M 0 

2.50 

2.44 — 2.61 

2.34 — 2.85 

^ 3:2 

-1 

-1 — -0.28 

-1 — 0.03 

^2:3 

0.11 

0.06 — 0.19 

-0.06 — 0.38 

Wl:l 

1 

0.39 — 1 

0.16 — 1 

case 3: 


a = 0.389 

a = 0.5 

a = 0.7 

M/Mg 

2.50 

2.44 — 2.61 

2.34 — 2.85 

^3:2 

-1 

-1 - -0.28 

-1 — 0.03 

^3:2 

-0.46 

-0.50 — -0.27 

-0.55 — 0.03 

Wl;l 

1 

0.39 — 1 

0.16 — 1 

case 4: 


a = 0.201 

a = 0.5 

a = 0.7 

M/Mg 

2.35 

2.19 — 2.61 

2.11 — 2.85 

^2:3 

-1 

-1 — 0.19 

-1 — 0.37 

^2:3 

-0.54 

-0.47 — 0.19 

-0.52 — 0.38 

<*>1:1 

1 

0.01 — 1 

0.10 — 1 


Table 2 Restriction on the parameters of the neutron star in the XTE J1701-407 source 
implied by the string loop oscillation model of HF QPOs. Presented values correspond to 
the hatched regions from Fig. 4. 


We assume relevance of resonance phenomena, e.g., a parametric resonance 
[41], at all of the three HF QPO events observed in the XTE J1701-407 source. 
We consider the rational frequency ratios vg : v T = 3 : 2 or vg : v T = 2 : 3, 
and vg ■ v T = 1 : 1, to be directly related to the observed values of the QPO 
frequencies in the XTE J1701-407 source. For all of the three observational 
events, we identify directly the frequencies vjj, v l with vg, v T or v T1 vg frequen¬ 
cies. There are four possible combinations of this identification enabled by the 
properties of the string loop oscillation model and presented in Table 1. In 
the cases 1, 2, all three of the resonant radii of the string loop oscillations, 
^ 3 : 2 , 7*1:1, t* 2 : 3 : are considered, while in the cases 3,4, the resonant radius ru, 
and two resonant radii r 3:2 (case 3), or two resonant radii r 2: 3 (case 4), are 
considered. 

In a given spacetime (fixed parameters M,a), the string loop oscillation 
model predicts for a string loop with a fixed parameter to that the resonance 
radii and the radial (or vertical) frequencies related to the upper observed 
frequency have to satisfy the relations 


r 3:2 < 11:1 < 7*2:3, ^3:2 > ^2:3 > ^1:1 


(23) 
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Fig. 5 Examples of the radial profiles of the string-loop oscillatory frequencies z/ r (thin 
curves) and vq (thick curves) as related to the three observational events are given for two 
representative situations allowed by the combination of the resonant points in the case 4. 
The parameters of the Kerr metric are mass M = 2.35Mq and spin a = 0.22 for the first 
row, M = 2.25Mq and spin a = 0.4 for the second row. The related values of the parameter 
lj are depicted in all the subfigures. Relevant resonant frequencies = 1153 Hz and 
^ 2:3 = (1110,740) Hz are given. The resonant radii 7*2:3 are close to the maximum of the 
radial profile of the radial mode of the string loop oscillations. 


where we denote as ^ m:n the upper frequency in the m : n resonance point - 
see Fig. 1. However, in the XTE J1701-407 source, there is (eqs. (1-3)) 


fl:l > /3:2, (24) 

therefore, it is impossible to explain all the observed resonances by assuming 
only one string loop with a fixed parameter to. We thus assume that at the 
resonant radii corresponding to the three observational events with HF QPOs, 
different string loops occur, being governed by different values of the parameter 
oo. The string loop parameter oo is thus considered to be a free parameter in a 
given spacetime with fixed parameters M, a. 

The procedure of fitting the string loop oscillation frequencies to the ob¬ 
served frequencies in the three observational events is presented in Figure 4 for 
all the four cases of possible combinations of the resonant radii of the string 
loop oscillations. At each of the three observed events, and each of the resonant 
points, the fitting is related to the upper of the observed frequencies (or the 
common frequency at rn); the precision of their measurement is also taken 
into account. The fitting procedure then gives for each of the observed events 
an allowed region of the parameter space of the spacetime parameters M, a, de- 
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u 


Fig. 6 Ratio of frequencies v T : vg calculated at the local maximum of i/ r (r) function 
determining the frequency of the radial oscillations of the string loops for different values of 
parameter a. 


termined by the limiting values of the string loop parameter oj £ (—1,1). Due 
to the degeneracy of the radial profiles of the string loop oscillation frequencies 
in the Schwarzschild spacetimes (a = 0), i.e., their independence of the stringy 
parameter w, the fitting predicts only one value of the mass parameter M for 
the spin a = 0 at each observational event. Extension of the allowed region 
related to the whole interval of string loop parameter oj £ (—1,1) (i.e., the 
interval of allowed values of M) increases with increasing spin a. Therefore, 
the string loop oscillation model implies a ’’triangular” limit on the spacetime 
parameters M, a for each of the observed events of HF QPOs - see Figure 4. 
In a given spacetime, the limits have to be satisfied simultaneously and we 
directly obtain the region of allowed values of the spacetime parameters, if the 
theoretical restrictions on the mass (M < 2.5Mq) and spin (a < 0.4 related to 
the Hartle-Thorne model of neutron stars) are also taken into account. (For 
completeness, we demonstrate in Figure 4 also the limits implied by the gen¬ 
eral restriction on the neutron star spin, a < 0.7). Along with the restrictions 
on the spacetime parameters M and a, we obtain simultaneously restrictions 
on the stringy parameter oj. Moreover, the restrictions on the radii where the 
resonant oscillations occur have to be also taken into the account, if we test 
the Hartle-Thorne models of the neutron star for concrete equations of state. 
The results representing the limits on the spacetime and string loop parame¬ 
ters M, oj are presented in Table 2 for the characteristic limiting values of the 
spin parameter a. 

Our results indicate that only the case 4, with both the resonant radii 
related to 7’ 2: 3, is physically realistic, as the other three cases, where the radius 
V 3-.2 enters the play, give unrealistic values of the mass and spin of the neutron 
star, and, moreover, the radius ?' 3,2 could be too close to the surface of the 
neutron star. Using the results obtained in the case 4, we can see that for the 
lowest value of the spin a = 0.2, the unique value of the neutron star mass 
reads M = 2.35M 0 . For the allowed spin a = 0.3, the mass parameter can be 
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in the interval M ~ (2.3 — 2.4)Mq, while for the upper limit on the spin of the 
source a = 0.4, the mass parameter can be in the interval M ~ (2.2 — 2.5)M 0 . 
Such large predicted values of the neutron star mass are in agreement with 
the assumption of applicability of the Kerr geometry in our study. Further, we 
obtained also restrictions on the stringy parameter to that has to be positive 
at the resonant radii rm, but negative at the resonant radii T 2 -. 3 - 

In Figure 5 we represent the situation of fitting the data from the three 
observational events in the case 4 for some realistic values of the spacetime 
and stringy parameters from the allowed interval of the parameters. Notice 
that the resonant points with the frequency ratio 2 : 3 are located in vicinity 
of the local maximum of the radial profile of the radial frequency. It means 
that the upper frequency can be nearly the same for a wide range of radii of 
the stable equilibria fulfilling the resonance condition at r 2 : 3 - The resonance 
can be thus obtained by slight shift in the radius, since the upper (radial) 
frequency remains nearly constant at vicinity of the local maximum of the 
radial frequency radial profile, while the vertical frequency varies strongly in 
this region. The situation is illustrated by Figure 6 demonstrating that the 
frequency ratio v T : ug is relatively close to 3 : 2 at the local maximum of ry (r) 
for negative values of the stringy parameter ui, but the ratio decreases for 
positive values of ui. The string loop oscillation model predicts negative values 
of the parameter u> at the resonance radii r 2:3 implying thus their occurrence 
in the vicinity of the local maximum of iy (r). 


6 Conclusions 

We have tested the string loop oscillation model introduced in [64] by fitting 
the frequencies of the three HF QPOs observed in the XTE J1701-407 LMXB 
source containing a neutron star, using the analytic formulae for the radial 
and vertical oscillations in the Kerr spacetime approximating the external 
geometry of the neutron star. This model, reflecting oscillations of the string 
loops governed by interplay of tension and angular momentum, gives relevant 
restrictions on the spacetime parameters M, a and the string loop parameter to 
that must be varied for the three observational events. The restricting results 
can be summarized as follows: 

— we cannot fit the observed data assuming only one string loop having a 
fixed value of the parameter ui, 

— the neutron star has to be rather fast rotating, having the dimensionless 
spin limited to the interval 0.2 < a < 0.4, 

— the neutron star has to be very massive, with mass parameter limited to 
the interval 2.1M 0 < M < 2.5M 0 . 

Since the XTE ,11701-407 neutron star has to be very massive, we can 
conclude that the application of the Kerr geometry in the fitting procedure 
is justified, as for the near-maximum-mass neutron stars the exterior Hartle- 
Thorne geometry has to be close to the exterior Kerr geometry, giving close 
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predictions of the physical phenomena occurring in their vicinity. Of course, 
we have to test in future, if application of the string loop model related to 
the Hartle-Thorne geometry with large values of q/a 2 > 3, when complex 
behaviour of the frequency relations is expected [79], could allow for fitting 
the data by neutron star models with intermediate or low mass, M < 2M 0 . 

In a future study, we would like to test the string loop oscillation model, 
potentially reflecting tension of magnetic field lines in the oscillating objects, 
also for the HF QPOs observed in typical atoll sources demonstrating large 
scatter of the observed twin HF QPOs. 

We can conclude that the string loop oscillation model of the HF QPOs 
occurring in the strong gravitational field of compact objects can fit both the 
high, kHz frequencies observed in the microquasars [64], i.e., binary systems 
containing a black hole, and the kHz frequencies observed in the special, low- 
luminosity atoll binary system XTE J1701-407 containing a neutron star and 
demonstrating a special set of HF QPOs untypical for the atoll sources. There¬ 
fore, it is worth to investigate the string loop oscillation model in more detailed 
way. First, we have to concentrate on the optical phenomena related to the 
radiation of the oscillating string loops that could be treated in the standard 
way introduced for the electromagnetic (gravitational) radiation of strings [84] 
and modification of the emitted radiation by strong gravity of a black hole or 
a neutron star. Second, we have to study the conditions for creation of the 
’’magnetic” string loops along the lines proposed in [21,22], giving the rela¬ 
tions of the physical parameters of magnetized plasma to the stringy angular 
momentum parameters J and ui. 
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